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It is shown numerically and analytically that the front propagation process in the framework of 
the Swift-Hohenberg model is determined by periodic nucleation events triggered by the explosive 
growth of the locahzed zero- eigenvalue mode of the corresponding linear problem. We derive the 
evolution equation for this mode using asymptotic analysis, and evaluate the time interval between 
nucleation events, and hence the front speed. In the presence of noise, we derive the velocity 
exponent of "thermally activated" front propagation (creep) beyond the pinning threshold. 
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One of fundamental models of pattern formation far 
from equilibrium is the Swift-Hohcnberg (SH) equation 
[Q. Despite relative simplicity of this model, and a se- 
rious limitation related to the fact that it is of the gra- 
dient type, it gives rise to a remarkably large variety of 
solutions. The SH equation has been intensively stud- 
ied in the past as a paradigm model of pattern forma- 
tion in large aspect ratio systems. [0J^]. More recent 
computations brought attention to the propagation of 
fronts between uniform stationary states of this equation 
and coarsening ||^. The computations have also demon- 
strated formation of stationary solitons, i.e. stable lo- 
calized objects in the form of a domain of one phase 
sandwiched inside another phase The interest was 
supported by applications of the SH model to marginally 
unstable optical parametric oscillators (OPO) 0|, and 
more numerical computations, leading to similar conclu- 
sions, were carried out in this context The SH model 
has also served as a convenient testing tool for the prob- 
lem of pattern propagation into an unstable trivial state 
[^,0. It is the purpose of this Letter to elucidate an- 
other aspect of pattern formation that can be modeled 
by the SH equation: coexistence between a pattern and 
a uniform state and stick-and-slip motion of interphase 
boundary which can be thought of as a particular kind 
of crystallization or melting. 

We shall write the basic equation in the form 
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At e > the trivial state u = undergoes supercriti- 
cal stationary bifurcation leading to a small-amplitude 
pattern with unit wavenumber k. The band of unstable 
wavenumbers widens with growing e, until it reaches the 
limiting value fc = 0, which signals the appearance of 
a pair of nontrivial uniform states, u = ±Ve ~ 1- The 
two symmetric states are stable to infinitesimal pertur- 
bations at e > |. At still higher values of e, a variety of 
metastable states become possible: (1) a kink separating 



the two symmetric nontrivial uniform states (Fig. ^); 
(2) a semi-infinite pattern, coexisting with either of the 
two nontrivial uniform states; (3) a finite patterned inclu- 
sion, sandwiched, either symmetrically or antisymmetri- 
cally, between semi-infinite domains occupied by nontriv- 
ial uniform states Q ; (4) a soliton (Fig. [^b) in the form of 
an island where one of the uniform states is approached, 
immersed in the infinite region occupied by the alterna- 
tive state; the basic soliton "atom" can be viewed as a 
single period of the pattern immersed in a uniform state. 
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FIG. 1. Stationary solutions to Eq. (|ll) for e = 1.78. (a) 
kink solution connecting u = ±y/e — 1; (b) ID localized so- 
lution (solid line) and radially-symmetric solution to Eq. (|^) 
(dashed line). 

The energy of the uniform state is higher than that 
of the regular pattern with the optimal wavenumber at 
e < 6.287 |l^]. In spite of the difference between the en- 
ergies of the uniform and patterned state, the interface 
between them remains immobile at moderate and large 
values e, and a multiplicity of localized states is linearly 
stable. The pattern propagation into a metastable uni- 
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form state or the reverse "melting" process at e ^ 1 is 
impeded by the self-induced pinning attributed to the os- 
cillatory character of the asymptotic perturbations of the 
uniform state, characterized by the complex wavenumber 
= —lzLiy/2e — 3. As e decreases, the stationary local- 
ized solutions lose stability and give rise to propagating 
solutions. 

We shall show in this Letter that the front propaga- 
tion process can be described in terms of periodic nu- 
cleation events triggered by explosive growth of the lo- 
calized zero-eigenvalue mode of the corresponding linear 
problem. Using asymptotic analysis, we shall derive the 
evolution equation for this mode which allows us to esti- 
mate the time interval between the nucleation events and 
hence the front speed. We shall also study the "thermally 
activated" front propagation (creep) at e > ec and derive 
the creep velocity exponent. 

We performed numerical study of Eq.(^ in one spa- 
tial dimension at various values of e. We found that 
the depinning transition for the single soliton occurs at 
e 1.74. The depinning threshold increases with the size 
of the patterned cluster, rapidly converging the limiting 
value e = Ec ~ 1.7574... for the semi-infinite pattern. Fig- 
ure H shows propagation of the semi-infinite pattern into 
a uniform stable state at e = 1.75. The front propaga- 
tion takes the form of well separated in time periodic 
nucleation events of new "atoms" of the "crystalline" 
state at the front. Between successive nucleation events, 
the solution remains close to the stationary semi-infinite 
pattern found at ec. This process resembles crystalliza- 
tion in equilibrium solids, with the important distinc- 
tion that the new "atoms" are created directly from the 
metastable "vacuum" state. The time between consecu- 
tive nucleation events diverges as e approaches the pin- 
ning threshold Figure ^ presents the average front 
speed as a function of e — Ec- This function can be fitted 
hyV = Vo^/ec ~ e, with Vq = 2.292. 
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FIG. 3. The average front speed as a function of the 
control parameter e near the critical point — 1.7574... 
Points represent numerics, the dashed line is the best fit 
V = 2.29(ec — e)^^^ and the dot-dashed line is the theoret- 



ical prediction V = 1.80(ec — e 
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At S — e ~ ec > 0, the immobile front solution uo{x, e) 
is linearly stable. Numerical stability analysis shows that 
there exists an exponentially decaying mode ui{x,e) lo- 
calized at the front with negative eigenvalue As which 
approaches zero as S-^^^. In addition, there is also an 
unstable front solution and a corresponding mode with 
positive eigenvalue A„ which also approaches zero as 6^^'^ . 
At the pinning threshold, these two solutions collide and 
disappear via a saddle-node bifurcation. The stationary 
front solution Uq{x) = uo{x,ec) and the corresponding 
mode structure Ui{x) = ui{x,ec) are shown in FigJ^. 
In the inset, we show the stable and unstable eigenval- 
ues Xs,u as functions of e. At e < ec, the front solution 
becomes non-stationary. Nonetheless, at \6\ <C 1, the nu- 
merical study suggests that the solution remains close to 
the stationary front solution Uoix) all the time except 
short intervals when a new roll nucleates. 




FIG. 2. The space-time plot of the switching wave prop- 
agation into the stable uniform phase ai e — 1.757. The 
horizontal axis is spatial coordinate, and the vertical axis is 
time. The system size is 200 (4096 mesh points), the time 
span is 200. 
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FIG. 4. The stationary structure Uo{x) (solid line) and 
the corresponding localized zero mode Ui{x) (dashed line) at 
e = Ec- Inset: eigenvalues of the stable and unstable localized 
modes of the stationary front solutions at e > tc- 

The process of front propagation can be analyzed 
within the framework of perturbation theory near the 
critical value e = ec- At \S\ 1, we can present the front 
solution in the form 

u{x,t) = Uoix) + \S\^^'^ui{x,t), 

where Uo{x) is the stationary front solution at S = 
and \S\^^^ui{x,t) is a small correction. This solution is 
uniformly valid at small positive S; it is also valid during 
quasi-stationary phases (away from nucleation events) for 
small negative S. Plugging this ansatz in Eq.(|^), we ob- 
tain 

dtUi = L[Uo]ui - |<5r/' [Woul - sign{5)Uo] + (2) 

where L[Uo] = Cc — 3C/q — (1 -I- V^)^ is the linearized 
SH operator at e^- In the lowest order in 5, Eq.(|^) 
yields the linear equation dtUi = L[Uo]ui. This equa- 
tion is always satisfied by the stationary translational 
mode U'q(x). In addition, the linearized operator L\Uq\ 
has a localized neutral eigenmode f/i (found numerically 
above, see Figure ^). Since all other eigenmodes have 
negative eigenvalues, the evolution of the system close to 
the bifurcation point can be reduced to single-mode dy- 
namics pd| ]. Therefore in the lowest order we can write 
ui{x,t) = a{t)Ui{x), where a{t) is the amplitude of the 
zero mode. 

Close to the bifurcation threshold dtUi = aUi can be 
treated as perturbation. Therefore, in the second order 
we derive 

L[Uo]u2 = aUi + |(5|i/2 [3a2{/o[/2 _ sign(5)[/o] . (3) 

Eq. (^ has a bounded solution if its r.h.s. is orthogonal 
to the zero mode Ui of the operator L. This results in 
the solvability condition for the amplitude a: 

ad=|J|i/2[sign(J)/3-7a2], (4) 

where 

/oo noo roc 

Ufdx, (3^1 UiUodx, 7 = 3 / UfUodx. 

At (5 > 0, a{t) reaches a stationary amplitude oq = 
{P/lY^'^- This value corresponds to the difference be- 
tween the stable front solutions at e and tc- At small 
negative 5, Eq.(^ describes explosive growth of a, which 
passes from -co to oo in a finite time = ira/ {\5\j3^)^^'^ . 
This explosion time gives an upper bound for the period 
between the nucleation events, after which the whole pro- 
cess repeats. The front speed is found as V = A/te, 
where A is the asymptotic spatial period of the pattern 



selected by the process of roll nucleation, and Te is the 
time interval between nucleation events. Our calculations 
give the value: V = 1.8|(5|^/^. This scaling is in a good 
qualitative agreement with the results of the numerical 
simulations, see Fig.^. However, the prefactor 1.8 is no- 
ticeably lower (about 25%) than the corresponding value 
Vb = 2.29 obtained by numerical simulation of Eq. (|l|). 

Let us discuss a possible reason for the discrepancy. 
Our numerical simulations show that the nucleation 
events produce slowly-decaying distortions behind the 
moving front. These distortions may effectively "pro- 
voke" a consequent nucleation event by creating an ini- 
tial perturbation of the zero mode Ui{x). It will lead 
to an increase of the front velocity. Although we have 
evidence for the importance of this effect, a systematic 
treatment of this process is very complicated and goes 
beyond perturbation theory. 

Effect of noise. Let us consider the effect of weak ad- 
ditive noise ri{x,t) in the r.h.s. of Eq. ([l]). For simplicity 
we assume that rj is delta-correlated noise with the in- 
tensity (temperature) T: 

{ri{x, t)ri{x', t')) = Td{x - x')S{t ~ t'). 

In the presence of noise, there is no sharp threshold 
for the onset of motion at e < Cc- Instead, thermally- 
activated motion (creep) occurs at e > Cc (see Fig. H). 
In this case the average creep velocity is determined by 
the intensity of noise. For e < Ec the noise will slightly 
increase the speed of the front. In contrast to the deter- 
ministic motion, the intervals between consecutive nucle- 
ation events are random. 




FIG. 5. The space-time plot of the switching wave prop- 
agation triggered by white noise with the temperature 
T = 0.00042 at e = 1.758 > Ec. The horizontal axis is spatial 
coordinate, and the vertical axis is time. The system size is 
200 (4096 mesh points), the time span is 1000. 

In order to estimate the effect of the noise at (5 > 0, we 
will treat rj as small perturbation. In this analysis, we 
shall not introduce scaling explicitly, since the scaling of 
noise that is comparable by its effect with deterministic 
perturbations cannot be determined a priori. Following 
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the lines of the above analysis, we project noise at the 
zero mode to obtain the solvability condition 



a a 



-1/2^; 



(5) 



where fj{t) = J^^'q{x,t)Ui{x)dx. For T = 0, Eq. (|) 
has a stable fixed point = \fWPl an unstable one 
o^u = ~\fWhl- At a < a^i, one has explosive growth of 
the solution to Eq. (^), while a — > at a > a„. Thus, we 
have to estimate the probability P for the amplitude of 
the zero mode a to be smaller than a„ . This quantity can 
be derived from the corresponding Fokker-Planck equa- 
tion for the probability density p(a, t) [T^ : 

.. = -^|^[(/^-7a>]+^5.^., (6) 

where we used (77^) = aT . This is the standard Kramers 
problem. The stationary probability p(a) is given by 



p exp 



2(53/2(/3a-7a3/3) 



T 



(7) 



The probability P is given by the integral P = 

we can use the saddle-point 
method, which gives the following result: 



P{a < au) 
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Since the time between the nucleation events t„ ^ ^/P, 
we find that the velocity of the front in the stable region 
is given by w l/r„ exp [-0.576^/^ /T] . 

At very large e > eo = 6.287, the flat state u = 
zby'e — 1 has lower energy than the periodic state. Nev- 
ertheless, the uniform state does invade periodic state at 
any e because of the strong self-induced pinning. With 
large-amplitude noise there will be some probability for 
the flat state to propagate towards periodic state by 
thermally-activated annihilation events at the edge of the 
periodic pattern, however t large e the probability of an- 
nihilation at the edge is of the same order as that in the 
bulk of the patterned state. Thus, for very large e and 
large T we may expect melting of the periodic structure 
both on the edge and in the bulk. 

The above results can be trivially extended to regu- 
lar two-dimensional periodic structures (rolls) selected by 
the SH equation . More interesting is the behavior of a 
2D hexagonal lattice which was predicted by weakly non- 
linear analysis near the bifurcation of nontrivial uniform 
solutions at e = I and found numerically . We antici- 
pate that propagation of the hexagonal structure into the 
uniform state will exhibit the same features of stick-and- 
slip motion as described above, and can be studied by 
similar methods. For SH equation (jl|) at e > |, hexag- 
onal lattices coexist with roll patterns which have lower 



energy, so the front propagation may in fact give rise to 
rolls. However, the hexagons may become dominant in 
the modified SH equation with an added quadratic non- 
linearity near e = 0, where the robust "crystallization" 
of the hexagonal lattice is expected. The work on this 
subject is now in progress. 
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